Symmetry fingerprints of a benzene single-electron transistor 
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The interplay between Coulomb interaction and orbital symmetry produces specific transport 
characteristics in molecular single electron transistors (SET) that can be considered as the finger- 
prints of the contacted molecule. Specifically we predict, for a benzene SET, selective conductance 
suppression and the appearance of negative differential conductance when changing the contacts 
from para to meta configuration. Both effects originate from destructive interference in transport 
involving states with orbital degeneracy. 



PACS numbers: 85.65.+h, 85.85.+j, 73.63.b 

Understanding the conduction characteristics through 
single molecules is one of the crucial issues in molecular 
electronics The dynamics of the electron transfer to 
and from the molecule depends on the intrinsic electronic 
spectrum of the molecule as well as on the electronic 
coupling of the molecule to its surroundings. 

In recent years the measurement of stability diagrams 
of single electron transistor (SET) devices has become a 
very powerful tool to do spectroscopy of small conduct- 
ing systems via transport experiments. Thus the capa- 
bility to perform three terminal measurements on single 
molecules 0, i, 0, 0, S, 0, E3] has been a fundamen- 
tal achievement for molecular electronics. Such molecu- 
lar transistors might display transport properties which 
are very different from those of conventional SETs. In 
fact, vibrational or torsional modes @, E3] and intrinsic 
symmetries/asymmetries of the molecule can hinder or 
favor transport through the SET, visible e.g. in the ab- 
sence/presence of specific excitation lines in the stability 
diagram of the molecular SET, or in negative differential 
conductance features. Many-body phenomena as e.g. the 
Kondo effect, have been observed as well 0, 00, E3]- 

Despite the experimental progress, the theoretical un- 
derstanding of the properties of single organic molecules 
coupled to electrodes is far from being satisfactory. On 
the one hand, numerical approaches to transport based 
e.g. on the combination of Green's function methods 
with density functional theory have become a standard 
approach to study transport at the nanoscale [l[. How- 
ever, this technique is not appropriate for the description 
of transport through a molecule weakly coupled to leads, 
due to the crucial role played by the Coulomb interaction 
in these systems. Hence, in |1 lj . an electronic structure 
calculation for a benzene molecule was performed in or- 
der to arrive at an effective interacting Hamiltonian for 
the 7r orbitals, to be solved to determine the I-V char- 
acteristics of a benzene junction. In this letter we con- 
sider the electronic transport through a benzene SET. 
Similar to 11], in order to devise a semi-quantitative de- 
scription, we start from an interacting Hamiltonian of 
isolated benzene where only the localized p z orbitals are 
considered and the ions arc assumed to have the same 
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FIG. 1: (color online) The two different setups for the benzene 
SET considered in this letter. 



spatial symmetry as the relevant electrons. The Hamilto- 
nian for the isolated molecule possesses 4 6 = 4096 eigen- 
states, to be calculated numerically, and whose symme- 
tries can be established with the help of group theory. 
Large degeneracies of the electronic states occur. For 
example, while the six-particles ground state (A\ g sym- 
metry) is non-degenerate, there exist four seven-particle 
ground states due to spin and orbital [Ei u symmetry) 
degeneracy. When coupling the benzene SET to leads 
in the meta and para configurations, these orbitals sym- 
metries lead to very different stability diagrams for the 
two configurations (see Fig.[2|). Striking are the selective 
reduction of conductance (Fig. [3|) and the occurrence of 
negative differential conductance (NDC) features when 
changing from para to meta-configurations. As shown 
in Fig. [4j the NDC effect occurs due to the formation 
of a blocking state at certain values of the bias voltage. 
The blocking is clearly visible by monitoring the position- 
dependent many-body transition probabilities which, at 
given values of the bias voltage, can exhibit nodes at the 
same position as one of the contacts. NDC for benzene 
junctions has been predicted also in [ll|, but in the para 
configuration and in presence of an external electromag- 
netic field. In our work NDC occurs despite the absence 
of an external field. Both the effects we predict originate 
from bias dependent interference of orbitally degenerate 
coherences, neglected in 



states: coherences, neglected in [111], are essential to cap- 
ture interference effects when solving the equations for 
the benzene's occupation probabilities. Interference phe- 
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nomena in transport through benzene have been recently 
discussed also in [12|, |l3[ . The parameter regime is how- 
ever very different, as both discuss the strong tunneling 
limit, where Coulomb blockade effects are not relevant. 

We start from the total Hamiltonian H = i?bcn + 
H loads + Ht, where the Hamiltonian for benzene reads: 

H hcn = & X 4a d ia + 6 X ( d L d t+la + 4+la d ta) 

+U X (ni T - 5) (n a - |) (1 ) 

i 

+V X (n iT + n 4 - 1) (n i+ i T + ra i+u - 1) . 

Here d\ a creates an electron of spin a in the p z orbital of 
carbon i, i — 1, . . . , 6 runs over the six carbon atoms of 
benzene and ni a = d\ a d i(y . This Hamiltonian is respect- 
ing the Deh symmetry of benzene and also the particle- 
hole symmetry. Mechanical oscillations are at this level 
neglected and all the atoms are considered in their equi- 
librium position. The parameters b, U, and V for isolated 
benzene are given in the literature [14| and are chosen 
to fit excitation spectra. Even if the presence of metallic 
electrodes is expected to cause a substantial renormaliza- 
tion of U and V, we do not expect the main results of this 
work to be affected by this change. The weak coupling 
suggests that the symmetry of the molecule will remain 
unchanged and with it the structure of the Hamiltonian 
([1]). The gate voltage V g is introduced by a renormal- 
ized on-site energy £ = £o — eVg and we conventionally 
set V g = at the charge neutrality point. We represent 
source and drain leads as two resevoirs of non- interact- 
ing electrons: ff lcads = J2akA e k ~ Pa)cl ka c akr7 , where 
a = L, R and the chemical potentials p a of the leads 
depend on the applied bias voltage Pl.r = Mo ± In 
the following we will measure the energy starting from 
the equilibrium chemical potential (1q = thus giving a 
negative energy to the holes in equilibrium. The coupling 



where T^r = ^\tL,fi^T^L y R equal the bare transfer rates 
with the constant densities of states of the leads T>l,r- 
Terms describing sequential tunnelling from and to the 



to source and drain leads is described by 

ft^Efic^+cUj, (2) 

aka 

where we define S aa as the creator of the electron in the 
benzene carbon atom which is closer to the lead a. In 
particular d Ra := d\ a ,d 5a respectively in the para and 
meta configurations, while d\ a := d\ a in both setups. 
Due to the weak coupling to the leads we can assume 
that the potential drop is all concentrated at the lead- 
molecule interface and is not affecting the molecule itself. 

Given the high degeneracy of the spectrum, the 
method of choice to treat the dynamics in the weak 
coupling is the Liouville equation method already used 
e.g. in [3, [l6[ . Starting point is the Liouville equa- 
tion for the reduced density operator a — Tri ea( i s {p| = 
— ^Trieads{[-ff, p]} where p is the density operator [171 ]. 
Due to the weak coupling to the leads we treat the ef- 
fects of i?x to the lowest non- vanishing order. The re- 
duced density operator a is defined on the Fock space 
of benzene but coherences between states with different 
particle number and different energy can be neglected, 
the former because decoupled from the dynamics of the 
populations, the latter being irrelevant due to their fast 
fluctuation (secular approximation). As a result we ar- 
rive at a generalized master equation (GME) where co- 
herences between degenerate states are retained. This 
approach is robust against the small asymmetries intro- 
duced in the molecule by the coupling to the leads or 
by deformation as far as the energy splitting that lifts 
the orbital degeneracy is comparable to the thermal en- 
ergy. The GME is conveniently expressed in terms of 
the reduced density operator a NE = Vne o'Vne, where 
Vne '■= Tli T \N E £t)(N E £t\ is the projection operator 
on the subspace of N particles and energy E. The sum 
runs over the orbital and spin quantum numbers £ and 
r, respectively. Eventually the GME reads 



(3) 



lead a are proportional to the Fermi function f(x—p a ) := 
fai x ) an d fa( x ) = 1 ~ fa ( x )> respectively. Still in the 
sequential tunnelling limit, but due to the presence of co- 



ne 



f+{H ben -E) + - Pa (H hen -E) 

7T 



dL* NE - 



dt 



fa (E — ffben) Pa(E — i?bcn) 

7T 



d aT a NE + h.c. 



E ^aV NE \dij+{E - E 1 ) a N - 1E 'd aT + d aT f-(E> - E) a N+1E ' S aT \v NE , 
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herences, also energy non-conserving terms are appearing 
in the generalized master equation, they are proportional 
to the function p a (x) = —Reip \ + 2tt( x ~ /•*«) where 
ip is the digamma function HiGj. Finally we write the 
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FIG. 2: (color online) Stability diagram for the benzene SET 
connected in the meta (above) and para (below) configura- 
tion. Red dot-dashed lines highlight the conductance cuts 
presented in Fig. [3] the red dashed line the region correspond- 
ing to the current trace presented in Fig. [4] The parameters 
used are U = 4|6|, V = 2.4|6|, T = 0.04|6|, T = 1(T 3 |6|. 

GME in the basis of the energy eigenstates for isolated 
benzene and find numerically the stationary solution. 

A closer analysis of the master equation allows also to 
define a current operator (one per molecule-lead contact) 



4=E VN E \d ar f+(H PPF - E)d) ar - 



NEt 



di T f-(E-H PPP )d aT V NE 



(4) 



and calculate the stationary current as the average 7l = 
Tr{cr sta t^L} = -Ir, with <7 sta t the stationary density op- 
erator. In Fig. [2] we present the stability diagram for 
the benzene SET contacted in the para (lower panel) 
and meta position (upper panel). Bright ground state 
transition lines delimit diamonds of zero differential con- 
ductance typical of the Coulomb blockade regime while 



a rich pattern of satellite lines represents the transitions 
between excited states. Though several differences can be 
noticed, most striking are the suppression of conductance 
and the appearance of NDC when passing from para to 
meta configuration. A zero bias cut of the stability dia- 
grams as a function of the gate voltage V g is plotted in 
Fig. [3] Only transitions between ground states are rel- 
evant for the conductance. The number of p z electrons 
on the molecule and the symmetry of the ground state 
corresponding to the conductance valleys are reported. 
The conductance in the meta and para configuration is 
the same for the N = 11 12 and N = 10 <-> 11 tran- 
sitions while it is systematically suppressed in all other 
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FIG. 3: (color online) Conductance of the benzene SET as 
a function of the gate voltage. Clearly visible are the peaks 
corresponding to the transitions between ground states with 
N and TV + 1 particles. In the low conductance valleys the 
state of the system has a definite number of particles and 
symmetry as indicated. Selective conductance suppression 
when changing from the meta to the para configuration is 
observed. 

cases. In other terms transitions between states with 
A or B symmetry, which do not have orbital degeneracy, 
are invariant under configuration change; transitions that 
involve an E symmetry, and hence imply twofold orbital 
degeneracy, are suppressed. Destructive interference be- 
tween orbitally degenerate states explains the systematic 
conductance suppression. By neglecting the energy non- 
conserving terms in ^ we derived an analytical formula 
for the conductance close to the resonance between N 
and N + 1 particle states: 



G 



N,N+1 



(AE) = 2e 5 



Li R 



EnmrW n \ d Lr\ N + 1 ' m )( N + m \ d Rr\ N ' n ) 



rL + r^ 



Y, nmaT \{X,n\d aT \N + l,m)\* 
I 



/'(AS) 



(S N+1 -S N )f(AE) + S N _ 



(5) 



where AE = E g: N~- E g ^+i + e V g is the energy difference 
between the benzene ground states with N and 7V+1 elec- 
trons diminished by a term linear in the side gate, n and 



m label the SAT-fold and ^Ar+i-fold degenerate ground 
states with N and N + 1 particles, respectively. Interfer- 
ence effects are contained in the numerator of the third 
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factor (overlap factor A), 
visible we remind that dL 



In order to make these more 
* Rr =^4A> where is the 
rotation operator of an angle <f> and <fi = n for the para 
while <f> — 27r/3 for the meta configuration. All eigen- 
states of -ffben are eigenstates of the discrete rotation op- 
erators with angles multiples of 7r/3 and the eigenvalues 
are phase factors. The overlap factor now reads: 



A = 



(N, n\d LT \N + l, m)|V*" 



(0) 



where <j) nm encloses the phase factors coming from the 
rotation of the states \N,n) and \N +l,m). Interference 
is possible only when Sn or > 1, that is in presence 

of degenerate states. It generates a considerable reduc- 
tion by passing from the para to the meta configuration 
as seen in Fig. [3l Interference also affects non- linear 




state is populated and the current is reduced. To visual- 
ize this, we introduce the probability (averaged over the 
z coordinate) 

P(x,y;t) = lim V) ± Az\(7glr\^ {r)\&g)\ 2 (7) 

L J-L/2 



for benzene to make a transition between the state \6g) 
and one of the states \7g£r) by adding or removing an 
electron in position f. Each of the lower panels of Fig. 
13] are surface plots of P{x, y; £) for the 7 particles ba- 
sis that diagonalizes the stationary density matrix at a 
fixed bias. The upper plot of the e panel describes the 
transitions to the blocking 7 particle state that accepts 
electrons from the source lead (close to the carbon 1) but 
cannot release electrons to the drain (close to carbon 5). 
The energy non-conserving rates prevent the complete ef- 
ficiency of the blocking by ensuring a slow depopulation 
of the blocking state. At large negative bias the blocking 
scenario is depicted in the panel a. We remark that only 
a description that retains coherences between the degen- 
erate 7 particle ground states correctly captures NDC at 
both positive and negative bias. 

To summarize, we analyzed the transport character- 
istics of a benzene based SET. The interplay between 
Coulomb interaction and orbital symmetry is manifested 
in a destructive interference involving orbitally degener- 
ate states, leading to selective conductance suppression 
and negative differential conductance when changing the 
contacts from para to meta configuration. 

We acknowledge financial support by the DFG within 
the research programs SPP 1243 and SFB 689. 



FIG. 4: (color online) Upper panel - Current through the 
benzene SET in the meta configuration calculated at bias and 
gate voltage conditions indicated by the dashed line of Fig. 
[21 A pronounced NDC is visible. Lower panels - Transition 
probabilities between the 6 particle and the 7 particle ground 
states for bias voltage values labelled a — e in the upper panel. 
The transitions to a blocking state is visible in the upper 
(lower) part of the e (a) panels. 

transport and produces in the meta configuration NDC 
at the border of the 6 particles state diamond (Fig. \2§ ■ 
The upper panel of Fig. [5] shows the current through 
the benzene SET contacted in the meta configuration 
as a function of the bias voltage. The current is given 
for parameters corresponding to the red dashed line of 
Fig. [2 In this region only the 6 and 7 particles ground 
states are populated. The 6 particle ground state is not 
degenerate. The 7 particle ground state is 4-fold degener- 
ate, though the twofold spin degeneracy is not important 
since spin coherences vanish in the stationary limit and 
the S z = 1/2 and —1/2 density matrices are equal for 
symmetry. At low bias the 6 particle state is mainly oc- 
cupied. As the bias is raised transitions 6 «-> 7 occur and 
current flows. Above a certain bias threshold a blocking 
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